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Abstract 


Cloud microphysics budgets in the tropical deep convective regime are analyzed based 
on a 2-D cloud resolving simulation. The model is forced by the large-scale vertical velocity 
and zonal wind and large-scale horizontal advections derived from TOGA COARE for a 20- 
day period. The role of cloud microphysics is first examined by analyzing mass-weighted 
mean heat budget and column- integrated moisture budget. Hourly budgets show that 
local changes of mass-weighted mean temperature and column-integrated moisture are 
mainly determined by the residuals between vertical thermal advection and latent heat 
of condensation and between vertical moisture advection and condensation respectively. 
Thus, atmospheric thermodynamics depends on how cloud microphysical processes are 
parameterized. 

Cloud microphysics budgets are then analyzed for raining conditions. For cloud- vapor 
exchange between cloud system and its embedded environment, rainfall and evaporation 
of raindrop are compensated by the condensation and deposition of supersaturated vapor. 
Inside the cloud system, the condensation of supersaturated vapor balances conversion 
from cloud water to raindrop, snow, and graupel through collection and accretion processes. 
The deposition of supersaturated vapor balances conversion from cloud ice to snow through 
conversion and riming processes. The conversion and riming of cloud ice and the accretion 
of cloud water balance conversion from snow to graupel through accretion process. Finally, 
the collection of cloud water and the melting of graupel increase raindrop to compensate 
the loss of raindrop due to rainfall and the evaporation of raindrop. 
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1. Introduction 


Cumulus parameterization has been one of major research issues for more than three 
decades in meteorological community since the atmospheric general circulation model 
(GCM) became a major research and operational tool. The main idea in the cumulus 
parameterization is to use the large-scale variables to estimate the amount of the precipi- 
tation and the net heating and moistening effects due to sub-scale disturbances (convection) 
that cannot be simulated in the GCM. Based on the observation that the tropical deep 
convection always occurs over the region of large-scale moisture convergence, Kuo (1965, 
1974) used the temperature difference between the cumulus cloud and the surrounding 
environment, and the large-scale moisture convergence to calculate the net convergence of 
moisture into the vertical column of the air of unit cross section produced by the large- 
scale circulation, and by evaporation from the earth surface. The large part of the net 
moisture convergence is assumed to produce precipitation. Based on the assumption that 
the rate of production of available potential energy by the large-scale processes is nearly 
balanced by the rate of consumption of the available potential energy by the convection, 
Arakawa and Schubert (1974) used the vertical distributions of the total vertical mass flux 
by the cumulus cloud ensemble, the total detrainment of mass from the ensemble, and 
the thermodynamical properties of detraining air to parameterize the cumulus convection. 
Such physically and observationally based cumulus convection parameterizations have been 
included in the numerical models which lead to success of simulations of atmospheric circu- 
lations. Different from the cumulus parameterization, cloud microphysics parameterization 
uses the cloud thermodynamic variables to directly predict the cloud variations based on 
the empirical relationships derived from laboratory experiments and in-situ observations. 
The cloud resolving models that include cloud microphysics parameterization schemes have 
demonstrated to simulate reasonable atmospheric thermodynamics in the tropics during 
the GARP (Global Atmospheric Research Program) Atlantic Tropical Experiment (GATE: 
e.g., Xu and Randall 1996; Grabowski et al. 1996) and Tropical Ocean Global Atmosphere 
Coupled Ocean- Atmosphere Response Experiment (TOGA COARE: e.g., Wu et al. 1998; 
Li et al. 1999). Willoughby et al. (1984), Lord et al. (1984), and Liu et al. (1997) showed 
that inclusion of ice microphysics in their numerical models led to more realistic simula- 
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tions of cloud structures of model hurricanes. To understand how the cloud microphysics 
parameterization improve simulations, it is necessary to know the cloud microphysics bud- 
gets for cloud water, raindrop, cloud ice, snow, and graupel, and understand the dominant 
conversion processes and their growth and decay mechanisms. 

In this study, the 2-D cloud resolving model is used to investigate the dominant 
cloud microphysical processes associated with tropical deep convection and their roles 
in atmospheric thermodynamics. In the next section, the cloud resolving model will be 
briefly described, and heat and moisture equations used in budget analyses will be derived. 
In section three, we will show that vertical potential temperature advection and latent 
heat of condensation, and vertical moisture advection and condensation are the largest 
terms, respectively, in the mass-weighted mean heat budget and in the column-integrated 
moisture budget. Further, the residuals between the two terms, respectively, accounts for 
hourly thermal and moisture variations, implying the fundamental importance of cloud 
microphysical processes in determining atmospheric thermodynamics. Thus, the cloud 
microphysics budgets are calculated to analyze the dominant cloud microphysical processes 
in the tropical deep convective regimes in section four. Scale analysis is conducted to 
explain dominant cloud microphysical processes. Summary is given in section six. 

2. Formulations for model, heat, and moisture budgets 


The cloud resolving model was originally developed by Soong and Ogura (1980), Soong 
and Tao (1980), and Tao and Simpson (1993). The 2-D version of the model used by Sui et 
al. (1998) and modified by Li et al. (1999) is used in this study. The governing equations 
with an anelastic approximation can be expressed as follows: 
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Here u, and w are zonal, and vertical air wind components; 9 and q v are air potential 
temperature and specific humidity respectively; C = (q c , qn 9ii qsi qg)i 9ci Qn Qi-> 5s? an( l 9g 
are the mixing ratios of cloud water (small cloud droplets), raindrop, cloud ice (small ice 
crystals), snow (density 0.1 g cm -3 ), and graupel (density 0.4 g cm 3 ), respectively; p is 
a mean air density which is a function of height only; wtv is a terminal velocity which is 
zero for cloud water and ice; n = {p/po) K K = ~ ' Ft is the gas constant, c p is the specific 
heat of dry air at constant pressure p, and p o =1000 mb ; c, e, d, and s denote condensation, 
evaporation, deposition, and sublimation, respectively; Q cn L v (c— e)+L s (d— s)+Lf(f— 

m) denotes the net latent heat release through phase changes among different cloud species, 
where / and m are fusion and melting, respectively; L v , L s , and Lf are heat coefficients 
due to phase changes; Qr is the radiative heating rate due to convergence of net flux 
of solar and infrared radiative fluxes; S c is source and sink of cloud content determined 
by cloud microphyical processes (see Appendix); Ds are dissipation terms; overbar ( ) 

5 


denotes a zonal-mean; subscript i, denotes an initial value, which does not vary with time; 
superscript ° denotes imposed observed variables in the model. 

To examine domain-mean heat and moisture budgets, the equations for T (— ird) and 
q v are derived from (6), 
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For mass- weighted mean heat budget, multiplying (7a) by p, and integrating the 
resulting equation vertically, and dividing it by the mass of the air column yield 
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and for column-integrated water vapor budget, multiplying (7b) by p, and integrating the 
resulting equation vertically yield 
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Here < () >= [Q]/[l], and [()] = J'* r p()dz, where z T is the height of top model level. H s , - 
P, and E are surface sensible heat flux, precipitation, and surface evaporation, respectively. 


3. Mean-weighted mean heat budget and column-integrated moisture budget 

The experiment analyzed in this study is conducted with the model forced by zonally 
uniform vertical velocity, zonal wind, and horizontal advections, which are derived by 
Professor Minghua Zhang who used the 6-hourly TOGA COARE observations within 
the Intensive Flux Array (IFA) region. The calculations are based on the constrained 
variational method on column-integrated budgets of mass, heat, moisture, and momentum 
proposed by Zhang and Lin (1997). Hourly sea surface temperature at the Improved 
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Meteorological (IMET) surface mooring buoy (1.75°S, 156 °E) (Weller and Anderson 1996) 
is also imposed in the model. The model is integrated from 0400 19 December 1992 to 0400 
9 January 1993 (local time). The horizontal domain is 768 km. A grid mesh of 1.5 km and 
a 12 second step are used in model integrations. More discussion of the model is reported 
in Li et al. (1999). Hourly simulation data are analyzed in the following discussions. 

Figure 1 shows the time evolution of vertical distribution of the large-scale atmospheric 
vertical velocity and zonal wind during 19 December 1992-9 January 1993 that are imposed 
in the model. Within 19-25 December 1992, upward motion was dominant, indicating- 
strong convection. From 26 December 1992-3 January 1993, downward motion became 
dominant, along with occasional upward motion, suggesting the dry phase. In the last 
few days, moderate upward motion occurred. The diurnal and two-day signals are also 
detected in Fig.la as indicated by Sui et al. (1997) and Takayabu et al. (1996) respectively. 
The large-scale westerly winds increase significantly in the lower- and mid- troposphere 
and reach their maximum of 20ms -1 at 600 mb around 3 January 1993 (Fig. lb). As 
mentioned previously, the model is also forced by the observed horizontal temperature and 
moisture advections (not shown), which have much smaller amplitudes than the vertical 
advections respectively. 

The comparison between simulation and observation is carried out by analyzing the 
linear correlation coefficients and the RMS of temperature and specific humidity, which aie 
shown in Fig. 2. The correlation can be as high as 0.7 for temperature only in the upper 
troposphere and near the surface, and it becomes lower in the mid- and lower-troposphere 
with the minimum of 0.1 around 730 mb. The correlation of specific humidity is about 0.7 
between 300 and 700 mb, but it reaches its minimum of 0.1 around 870 mb. The RMS 
of temperature increases from 0.7 °C near the surface to 2.5 °C around 200 mb, whereas 
the RMS of specific humidity is 0.15-0.75 gkg ~ l . Some explanations of differences between 
simulation and observation can be referred to Li et al. (1999). 

Fig. 3 shows the each term of mass- weighted mean heat budget and column-integrated 
moisture budget versus local changes of mean-weighted mean temperature and column- 
integrated moisture (precipitable water) respectively. In the heat budget, the local changes 
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of mass- weighted mean temperature are mainly caused by the sum of latent heat of conden- 
sation and vertical potential temperature advection, because of small horizontal thermal 
advect.ion and surface sensible heat flux, and nearly constant radiative cooling (Fig. 3a). In 
the moisture budget, the local changes of precipitable water are mainly due to sum of pre- 
cipitation and vertical moisture advection, since the small horizontal moisture advection, 
and nearly constant surface evaporation (Fig. 3b). 

Fig. 4 shows that the latent heat of condensation and vertical potential temperature 
advection in mass-weighted mean heat budget and the precipitation and vertical moisture 
advection in column-integrated moisture budget have same orders of magnitudes but the 
opposite signs, indicating these terms cancel in large parts. To examine the dependence of 
the residuals between these largest terms on intensity of convection, we plotted the sum 
of mass-weighted mean latent heat of condensation and vertical potential temperature 
advection, and sum of condensation and column-integrated vertical moisture advection 
versus rain rates in Fig. 5. The thermal sum increases from small negative values to 
relatively large positive values as the rain rate increases (Fig. 5a), whereas the moisture 
sum decreases from small positive values to relatively negative values as the rain rate 
increases (Fig. 5b). Atmosphere tends to be warmed and dryed when surface rain rate is 
larger than 30 mmday . 

Figs. 3-5 indicate that the cloud microphysics parameterization is very crucial to 
determine the atmospheric thermodynamics. The dominant effect of sum of latent heat 
of condensation and vertical potential temperature advection on temperature variation is 
also evident in the mid-troposphere (309, 487, and 694 mb) as shown in Fig. 6. In the 
upper-troposphere (194 mb), the increase of temperature is mainly due to vertical heat 
flux convergence, whereas the decrease of temperature results from radiative cooling (Fig. 
6a). In the lower-troposphere (907 mb), vertical heat flux convergence, the latent heat of 
condensation, and vertical potential temperature advection are responsible for temperature 
variations (Fig. 6e). 

The dominant effect of sum of condensation and vertical moisture advection on mois- 
ture variations is only seen at 487 mb (Fig. 7c). In the upper-troposphere (194 mb and 309 
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mb), the condensation and vertical moisture advection and vertical moisture flux conver- 
gence cancel each other (Figs. 7a and 7b). In the lower-troposphere (694 mb and 907 mb), 
vertical moisture flux convergence is an important factor for the local moisture changes 
(Figs. 7d and 7e). Although the sum of condensation and vertical moisture advection 
is the dominant factor in column- integrated moisture budget, the vertical moisture flux 
convergence is the important process at the individual vertical levels. 

4. Dominant cloud microphysical processes 


Before discussion of dominant cloud microphysical processes, column-integrated cloud 
contents versus surface rain rate and vertical structures of mean cloud contents are ex- 
amined (Fig. 8). Cloud water, raindrop, graupel, and total cloud contents increase with 
increasing surface rain rate, whereas cloud ice and snow are much less sensitive to surface 
rain rate. Maximum of total cloud contents is between 500 and 600 mb, which is largely 
contributed by graupel, raindrop, and cloud water. Maxima of cloud ice and snow appear 
around 250 mb and 350 mb levels respectively. 

There are two aspects of cloud microphysical processes: cloud- vapor exchange between 
the cloud system and its embedded environment, and conversions among cloud contents 
inside the cloud system. In this study, notations for cloud microphysics are consistent with 
the previous studies. The rate of conversion for cloud microphysics is always denoted by a 
positive value, and the sign before the conversion term represents the direction conversion 
goes to. For cloud- vapor exchange between cloud system and the environment, there are 
seven cloud microphysical processes (see Ala). Only three of them are important (Fig. 
9a). Two sinks for vapors are conversion from vapor to cloud water by the condensation of 
supersaturated vapor (- Pcnd ) and conversion from vapor to cloud ice by the deposition 
of supersaturated vapor (- Pdep )• One source for vapor are conversion from raindrop to 
vapor by evaporation of raindrop ( Prevp )• Vertical profiles of mean cloud microphysical 
processes (Fig. 9b) display that the maxima of Pcnd and Pdep are at 600 mb and 350 mb 
levels respectively. The Prevp increases from 500 mb to 700 mb, and keep the constant 
below 700 mb level. 
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Inside the cloud system, there are nine cloud microphysical processes for cloud water 
(see Alb). Three of them are dominant (Fig. 10). Two sinks for cloud water are growth 
of raindrop by the collection of cloud water (- Pracw ) and growth of graupel by the 
accretion of cloud water (- Pgacw )• One source for cloud water is growth of cloud water 
by the condensation of supersaturated vapor ( Pcnd )• Although it is relatively small, 
growth of snow by the accretion of cloud water (- Psacw ) is also included to close the 
budget. Since the sum of these cloud microphysical processes are much smaller than the 
dominant processes, the condensation of supersaturated vapor maintains conversion from 
cloud water to raindrop, snow, and graupel by the collection and accretion of cloud water 
respectively. The large positive values of Pcnd are between 550 mb and 850 mb, whereas 
the large negative values of Pracw and Pgacw are around 600-800 mb and 550 mb level 
respectively. Four terms cancel each other in large part (Fig. 15a). The sum of these 
terms show maximum positive values at 900 mb level and maximum negative value at 550 
mb, and their magnitudes are about the same so that column-integrated sum is neai zero. 

For raindrop, there are twelve cloud microphysical processes (see Ale). Three of 
them are important (Fig. 11). Two sources for raindrop are the collection of cloud wa- 
ter (Pracw) and the melting of graupel ( Pgmlt )• O nc sink for raindrop is growth of 
vapor by evaporation of raindrop (- Prevp ) The sum of the rates of these processes is 
linearly proportional to surface rain rate. Thus, rainfall and evaporation of raindrop are 
compensated by melting of graupel and collection of cloud water. The maximum of cloud 
microphysical processes related to raindrop is at about 570 mb level, which is mainly con- 
tributed by Pgmlt, as well as Pracw (Fig. 15b). Pracw and -Prevp have the similar 
vertical structures but have the opposite signs. Thus, Pracw and -Prevp cancel each 
other in large parts. 

For cloud ice, there are nine cloud microphysical processes (see Aid). Three of them 
are important (Fig. 12). Two sinks for cloud ice are growth of snow by the conversion 
of cloud ice (- Psaut ) and by the riming of cloud ice (- Psfi )• One source for cloud 
ice is the deposition of supersaturated vapor ( Pdep )■ Since the sum is about zero, the 
deposition of supersaturated vapor maintains conversion from cloud ice to snow through 
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conversion and riming of cloud ice. Maximum positive value of Pdep occurs at 350 mb 
level, whereas maximum negative values of Psaut and Psfi appear at 400 mb and 250 
mb levels respectively (Fig. 15c). Maximum positive and negative values of the sum are 
at 200 mb and 300 mb levels respectively, and they have similar magnitudes that lead to 
a near-zero column integration. 

For snow, there are fifteen cloud microphysical processes (see Ale). Four of them are 
important (Fig. 13). Three sources for snow are the conversion of cloud ice {Psaut), the 
riming of cloud ice (Psfi), and the accretion of cloud water ( Psacw , note 84 = 1 here). 
One sink for snow is growth of graupel by the accretion of snow (- Pgacs )■ Since the 
sum is almost zero, the growth of graupel by the accretion of snow is nearly balanced by 
the conversion from cloud ice to snow through the conversion and riming processes, and 
conversion from cloud water to snow through the accretion process. Maximum positive 
values of Psaut, Psfi, and Psacw are at 400 mb, 250 mb, and 500 mb levels respectively. 
Maximum negative value of -Pgacs is at 370 mb level with the second maximum negative 
value at 530 mb level (Fig. 15d). The sum has a maximum positive value at 250 mb level, 
which is mainly contributed by Psfi, and a maximum negative value at 570 mb, which is 
mainly contributed by - Pgacs ■ Column- integrated sum is about zero. 

For graupel, there are fourteen cloud microphysical processes (see Alf). Three of them 
are important (Fig. 14). Two sources for graupel are the accretion of cloud water {Pgacw ) 
and the accretion of snow {Pgacs)- O ne sink for graupel is growth of raindrop by melting 
of graupel (- Pgmlt )■ Since the sum is much smaller than the dominant processes, melting 
of graupel is compensated by the accretion of cloud water and snow. Vertical structures of 
mean cloud microphysical processes (Fig. 15e) show that the sum has a maximum negative 
value at 570 mb level, which is contributed by -Pgmlt, and maximum positive value at 
500 mb level, which is contributed by Pgacw and Pgacs- Thus, the column-integrated 
sum becomes very small. 

Figure 16 summarizes mean column-integrated cloud microphysical budget. For ex- 
change between cloud system and the atmospheric environment, rain rate (10.6 mm day 1 ) 
and evaporation of raindrop (6.4 mm day -1 ) are nearly compensated by the condensation 
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(14.0 mm day -1 ) and deposition (2.5 mm day- 1 ) of supersaturated vapor. Inside the cloud 
system, the condensation of supersaturated vapor (14.0 mm day- 1 ) supports conversion 
from cloud water to raindrop (9.4 mm day- 1 ), snow (0.9 mm day ~ l ), and graupel (3.3 mm 
day -1 ) through the collection and accretion processes. The deposition of supersaturated 
vapor (2.5 mm day- 1 ) supports the conversion from cloud ice to snow through the conver- 
sion (1.6 mm day -1 ) and riming (0.7 mm day- 1 ) processes. The conversion and riming of 
cloud ice and the accretion of cloud water maintain the conversion from snow to graupel 
by the accretion process (3.0 mm day _1 ). Finally, the collection of cloud water (9.4 mm 
day -1 ) and the melting of graupel (6.5 mm day -1 ) increase raindrop to balance the loss 
of raindrop due to precipitation (10.6 mm day 1 ) and evaporation of raindrop (6.4 mm 
day- 1 ). 


5. Scale analysis 


The source for development of clouds stems from the difference between atmospheric 
specific humidity and saturation specific humidity for water and ice. The schemes of Tao 
et al. (1989) for calculations of Pcnd (A7) and Pdep (A17) indicate that the sum of 
Pcnd and P DE p is linearly proportional to the difference and is modified by temperature 
and mixing ratios of cloud water and ice. Pcnd and Pdep are partitioned by the linear 
function of temperature. Reference temperatures are T 0 o for Pcnd and T 0 for Pdep so 
that magnitude of Pcnd is much larger than that of Pdep ■ Th e major sink of clouds is 
Prevp (A5). For the evaporation and deposition processes (A3, A4, A5, A24, A30), the 
second terms are larger than the first terms. Thus, the ratios of [Prevp] to [ Pm i.'j'c] • 
[Prevp] to [Pmlts], [Prevp] to [Psdep], and [Prevp] to [Pgdep] can be estimated by 
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In estimation of (9), the fact, that the covariance between 5 — 1 and q r is larger 
than the covariances between 5 — 1 (5» — 1) and q s and between 5—1 (Si — 1) and q g 
is considered. The raindrop has larger falling velocity than snow and graupel. Thus, 
(^-)2 ~ 12 . 5 , (^r)a ~ 51.3. Therefore, [Prevp] is more than one order of magnitude 
larger than [ Psdep ], [Pmlts], [Pgdep], an d [ Pmltg ]• 

For the melting processes, the second terms of Pgmlt (A8) and Psmlt (A9) are 
much larger than the first terms. The ratio of [Pgmlt] to [Psmlt] can be expressed by 

[Pgmlt] ,<i i[{T-Ta)(-f) 2 ] , ^ 

[wT^ [(r _!■„)(*)***]• ( } 


The graupel has larger falling speed and mixing ratio than does the snow. Thus, ( fl , ) 2 
4.1. The magnitude of q g is twice as large as that of q s . Therefore, [Pgmlt] is about one 
order of magnitude larger than [Psmlt]- 

For the accretion and collection processes, the second terms of Praci (A10), Pracw 
(All), and Pjacr (A15) are much larger than the other terms. The ratio of [Pracw] to 
[Praci] becomes 

[Pracw] ^ [p~*qcQr] 

[Praci] [p~*qiQr] 

q c and q r are large in the mid and lower troposphere, whereas q x are large in the upper 
troposphere. Thus, the covariance between q c and q r is much larger than the covariance 
between qi and q r . [Pracw] is at least one order of magnitude larger than [Praci] ■ 
Similar analysis shows [Pracw] is one order of magnitude larger than [Pi AC r] ■ For the 
accretion processes, Pgacs (A27) is proportional to the covariance between q s and (],/ . 
whereas P RAC s (All), Psacr (A23), and Pqacr (A26) are proportional to the covariances 
between q r and q s and between q r and q g . Since the covariance between q s and q g is much 
larger than the covariances between q r and q s and between q r and q g , [Pgacs] is much 
larger than other three processes. 
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For the accretion, collection, and riming processes, the collection efficiency coefficients 
of Psacw (A20), Pgacw (A28), and Pwacs (A29) are one order of magnitude larger than 
those of Psaci (A19) and Pqaci (A25). The falling speed coefficient of the graupel (a) 
is more than one order of magnitude larger than that of the snow (a"). Thus, [Pgacw] I s 
the largest process compared to other four processes. 

Since the second term is larger than the first term in Psaut (A18), and Ati has 
the similar magnitude as At in Psfi (A22), [Psaut] and [ Psfi ] have same orders of 
magnitudes. Psaut and Psfi occur when T < T 0 , whereas p IMLT (A16) occurs when 
T > T 0 . [-ATfai m“ 2 + npq c r 2 U ir )] of P S fw (A21) is very small (not shown). Therefore, 
magnitudes of [Psaut] and [Tsf/] are much larger than those of [ Pimlt ] and [ Psfw ]• 
Since cloud ice hardly exceeds the threshold, Praut (A13) is small. Mixing ratio of cloud 
water is small when T < Too- Thus, [Pihom] is negligible. Calculation indicates that 
|_2iHo_ e o.5|r— t 0 | | < o.Ol. [Pidw] is at least one order of magnitude smaller than [Psfw]- 

6. Summary 

Dominant cloud microphysical processes and their roles in atmospheric thermodynam- 
ics are investigated by using a 2-D cloud resolving model. The model is integrated for 20 
days under the forcing of large-scale vertical velocity and zonal wind, as well as the large- 
scale horizontal advections derived from TOGA COARE data. Analyses of the hourly 
mass-weighted mean heat budget and the hourly column-integrated moisture budget show 
that the local changes of temperature and moisture are mainly caused by residuals between 
vertical thermal advection and latent heat of condensation and between vertical moisture 
advection and condensation respectively. This indicates that the cloud microphysical pro- 
cesses play important roles in determining atmospheric thermodynamics. 

Two important aspects of cloud microphysical processes are cloud-environment inter- 
action and interactions among cloud contents. The cloud-environment interaction is rep- 
resented by cloud- vapor exchange between cloud system and its embedded environment. 
This exchange eventually affects the development of convection and atmospheric thermo- 
dynamics. The mean cloud microphysics budget shows that precipitation and evaporation 
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of raindrop are largely compensated by the condensation and deposition of supersaturated 
vapor. The interactions among cloud contents are represented by the conversion among 
cloud contents, which are associated with collection, conversion, accretion, and riming 
processes. The condensation of supersaturated vapor converts vapor into cloud water. 
The collection and accretion of cloud water further enhance raindrop, snow, and graupel 
respectively. The deposition of supersaturated vapor converts vapor into cloud ice. The 
conversion and riming of cloud ice and the accretion of cloud water generate snow, and 
the accretion of snow produces graupel. Finally, the collection of cloud water and melting 
of graupel balance the loss of raindrop due to rainfall and its evaporation. 
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Appendix 

Microphysical processes parameterized in cloud resolving model 


Microphysics parameterizations in the cloud resolving model used in this study are 
based on the schemes proposed by Rutledge and Hobbs (1983, 1984; referred to as RH83 
and RH84 thereafter), Lin et ah, (1983; LFO), Tao et al. (1989; TSM), and Krueger et al. 
(1995; KFLC) respectively. Corresponding equations are as follows. 

Sq v = -Pcnd - Pdep ~ (1 - tfi )Psdep{T < To) - (1 - 8i)P G dep(T < T 0 ) 

+Prevp + Pmltg{T > To) + Pmlts{T > To ), (Ala) 

S qc = -Psacw ~ Praut - Pracw - Psfw(T < T 0 ) - Pgacw + Pcnd 

—Pihom{T < Too) + Pimlt{T > To) — Pidw(Too < T < To), (^41^) 

Sq r = Psacw {T > To) + Praut + Pracw + Pgacw (T > To) — Prevp 

+Pracs(T > To) — Piacr(T < To) — Pgacr(T < To) — Psacr(T < To) 

-Pgfr(T < To) + Psmlt{T > To) + Pgmlt{T > To), (^41°) 

S, H = —Psaut{T < To) — Psaci(T < To) — Praci(T < To) — Psfi(T < To) 

-Pgaci{T < To) + Pihom(T < T 0 o) - Pimlt{T > T 0 ) + Pdep 

+Pidw(Tqo < T < To), (Aid) 

S qs = Psaut{T < T 0 ) + Psaci(T < T 0 ) + S a P S acw{T < T 0 ) + Psfw{T < To) 
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+Psfi(T < To) + SoPraci{T < To) — Pracs(T > To) — Pgacs 
-Psmlt{T > T 0 ) - (1 - S 2 )Pracs(T < To) + 8 2 Psacr(T < T 0 ) 

+(1 — Si)Psdep(T < To) - Pmlts(T > To) + <5 3 Piacr{T < T 0 ) 

— (1 — <$4 )Pwacs(T < To), (Ale) 


S qg - (1 - 5 3 )Praci(T < To) + Pgaci(T < T 0 ) + Pgacw(T < To) + Pgacs 

+(1 - ^4 )Psacw(T < T 0 ) + (1 - 83 )Piacr(T < To) + Pgacr(T < To) 
+Pgfr(T < To) + (1 — S 2 )Pracs(T < To) + (1 — 8^)Pw acs(T < To) 
—Pgmlt(T > To) + (1 — 8 i)Pgdep(T < T 0 ) — Pmltg(T > T 0 ) 

+(1 — 8 2 )Psacr(T < T 0 ), (Alf) 

where 


| 1 if q c + qi > 10 8 gg 1 ,T < T 0 , 
\ 0 otherwise, 


(A2a) 


82 


1 if q s + q r < 10 4 gg 1 ,T < T 0 , 
0 otherwise, 


(. A2b ) 


| 1 if q r > 10 4 gg 1 ,T < T 0 , 
\ 0 otherwise, 


= 


1 if q s < 10 ~ 4 gg-\q c > 5 x 10- 4 ^ _1 ,T < T 0 , 
0 otherwise, 


(A2d) 


T 0 = 0 °C, and T 00 = -35 °C. The microphysical processes in the terms of the right-hand 
side of (Al) and corresponding schemes are described in Table Al. 

The mathematical formula of cloud microphysical proccsess are shown as follow: 
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Pmltg = 


2itNog(S — 1) r 0.78 , f ap.i.po.i T(6/2 + 5/2), 

^ + bo ~ [ aT + 031( 7 ) ( 7> 

G 


(-43) 


where Nog{= 4 x 10 6 ra -4 ) is the intercept value in graupel size distribution; S is the satu- 
ration ratio with respect to water; a(= 19.3 m 1- *\s -1 ) is the constant in fallspeed relation 
for graupel; b(= 0.37) is the fallspeed exponent for graupel; A! = ~ 1); B' — 

— ryAL — (Pruppacher and Klett 1978); K a (= 2.43x 10 -2 Jm~ 1 s~ 1 K 4 ) is the thermal con- 
ductivity coefficient of air; M w = 18.0160 is the molecular weight of water; x(= 2.26 x 10 
m 2 s _1 ) is the diffusivity coefficient of water vapor in air; B.(= 8.314 x 10 3 Jkmol~ l K 4 ) 
is the universal gas constant; e ws (= Nm~ 2 ) is the saturation vapor pressure for water; 

[= ( IpgNog ) j j s f+e slope of graupel size distribution; Pg{ = 400 kgm~ 3 ) is the density 
of graupel; p( = 1.718 x 10 — ° kgm~^ s~ l ) is the dynamic viscosity of air; T is the Gamma 
function. 


Pmlts = 


4iVos(S-l) f 0.65 
p(A> + B ') 1 A 2 


+ 0.44( 


2 -£) 2(g»)i 

p p 


r(V2 + 5/2) 

\ b/2+5/2 J’ 
A S 




where N os {= 4 x 10 6 m~ 4 ) is the intercept value in snowflake size distribution; a"{= 1.139 
rn 1-b s -1 ) is the constant in fallspeed relation for snow; b(= 0.11) is the fallspeed exponent 
for snow; As[= ( URsBos ) ] is the slope of snow size distribution; ps{= 100 k:gm~ 3 ) is the 
density of snow. 


Prevp = 


27tNqr(S — 1) .0.78 
P(A' + B>) [ Xl 


,o!p, i /Po x 1 r(3) 

+ 0.31 (-P )H-)4^r], 

P P X R 


(^45) 


where N 0R (= 8 x 10 6 m -4 ) is the intercept value in raindrop size distribution; a'(= 3 x 10 3 
s _1 ) is the constant in linear fallspeed relation for rain; Xr[= is the slope of 

snow size distribution; Pl(— 10 3 kgm~ 3 ) is the density of snow. 


PlMLT 


Qi_ 

A V 


("46) 


where At is the time step. 
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PCND = 


1 T-Too 


q v - (q ws + qis ) 


AtTo-Tool + i^ + ^K^f^ + ^^y 

where q ws and q ts are the saturation mixing ratio for water and ice respectively; A\ 
( r- 7 35 .S^ ; A 2 = (j 6 - 5 7 5 6 > ; B\ = 17.2693882; B 2 = 21.8745584. 


m 


— 2tx 


r 0.78 


— , . r w,iu _ / a/9 1 , po x i r(fe/2 + 5/2) , 

^ + o-3i<f >'<7>‘ * - 1. 


M8) 


— 27 r r 0.65 , (i n p N i / /^o \ i r(fe /2 + 5 / 2 ) . 

Psw = — if.(T - + ».«(y )'(;) ' -7^ 


pL/ 


A| 


a: 


(A9) 


*r „ „ ,/>„ a ,a 0 r(3) a,r(4) a 2 r(5) a 3 r(6) 

PrACI = -jqiERlNoR{ — ) 2 [ — 73 1 74 I p r , 6 J, 

4 p a r a r a r a d 


(-410) 


'ft 


where E RI (= 1) is the rain/cloud ice collection efficiency; a Q = -0.267 ms -1 , « x = 
5.15 x 10 3 s“\ a 2 = -1.0225 x 10 6 a 3 = 7.55 x 10 7 m" 2 *" 1 , which are the 

coefficients in polynomial fallspeed relation for rain. 


7T 


Pracw = —qcE R cNo R { 

4 P 


po,ir(i 0 T(3) «i.r(4) a 2 £(5) a3£(6) 


A 3 + A 4 + A 5 + A 6 ^ (4H) 

a r a r a r a r 


where E R c(= 1) is the rain/cloud water collection efficiency. 


PrACS = EsrR 2 ~\Vr — V s\NqrNqs{^—) 2 [ * 6 \ 4 \5\2 u \3 1’ (412) 

p p a s Ar a s a r a s a r 

where E SR (= 1) is the snow/rain collection efficiency; Vfl[= (-0.267 + - 2 04 /^ 10 - + 

9.06xio 3 )(go)|] j g j-j ie mass- weighted fall-speed for rain; V s[= a" (^)^j the rnass- 

*R.P S' 

weighted fall-speed for snow. 


Praut = «( 7 c - tfo), 


(413) 
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where a(= 10 _3 s _1 ) is the rate coefficient for autoconversion; c/o(= 1-25 x 10 3 gg x ) is 
the mixing ratio threshold. 


Pidw — 


n 0 e 


0.5|T-T 0 | 


1000/9 


M 


(HP 


noe 0.5\T-T 0 \} 


l>2 


(AU) 


where n 0 = 10 -8 m -3 ; b 1 and b 2 are the positive temperature-dependent coefficients tabu- 
lated by Koenig (1971). 


Piacr = 


iERI ^^N 0R (^) 

1 24 p v p ’ 


r a 0 r(6) , ftl r(7) , a 2 r(8) o 3 r(9) . 

[ \ a i I »G I 




A 7 k 


^ R 


4 


(A15) 


where n r , (= £2?-) is the number concentration of cloud ice crystals; M t (= 6 x 10 12 kg) is 

clK Mi' 

the average mass of a cloud ice particle. 


PlHOM = 2^’ 


(7116) 


Pdep — 


1 T 0 -r 

A(T„-T 00 l + (^ 


( iv {Qws + Qis) 

i A 2 qtqis \( L v T-7qq 

' qc+qi ' ' Cp To — Too 


+ 


L s T 0 -T 
c p Tq — Too ' 


(7117) 


PsAUT 


P<li 


p At 


O.G(T-Tn) 


7 


where M max {= 9.4 x 10 ~ l0 kg) is the maximum allowed crystal mass. 


(7118) 


0 ira"qiEsiNos f P0\h r(6 + 3) 

P S ACI= J Ip"’ 

where Esi{= 0.1) is the snow/cloud ice collection efficiency. 


„ Ka"q c EscNas f Po x i T(b + 3) 

psacw- — 3 (7) ip-- 

where E S c{= 1) is the snow/cloud water collection efficiency. 


(7119) 


(t120) 
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PSFW ~ — 7T~ (^1 m ir + KpQc^Uir), 

JUirAti 


(A 21 ) 


where At x [= ^ - mj"* 3 )] is the timescale needed for a crystal to grow from 
radius r 0 to radius r; m ir {= 3.84 x 1Q~ 9 kg) and U ir {= lms" 1 ) are the mass and terminal 
velocity of an ice crystal r(= lOO^ra); m,i r ( = 2.46 x 10 w kg) is the mass of an ice crystal 
7'o(= 40 /tm). 


Psfi = 


Qi 

Ah' 


(A 22) 


PsACR = E S rTT 2 — \Vs - Vr\No R Nos {'-^) 2 [ , 6 . ' + 

P P A R A S 


P 0 


A 5 «A1 


+ 


0.5 1 

A/jA|' 


(A23) 


PSDEP = 


4iVos(5i - 1) 0.65 
p(A» + B ") [ A| 


+ 0.44( 


P P 


1 

4 


r(b/2 + 5/2) , 

. b/2+5/2 1’ 


where 5^ is the saturation ratio with respect to ice; A" = ( L ^ w - 1); -E" 

e is (= Nm~ 2 ) is the saturation vapor pressure for ice. 


(A24) 

RT . 
xM w e is 5 


Pgaci = 


TuigiEgiNpo 

4 



r(fr + 3 ) 
A^ 3 ’ 


where Egi{= 0.1) is the graupel/cloud ice collection efficiency. 


(2425) 


PgACR = EgrK 2 — \Vg - Vr\NorNqg ( ( -^) 2 [-rgr h \ 5 T 2 + \ 4 \ 3 ]’ (A 26 ) 

P P A R A G a r a g a r a g 

where Egr(= 1) is the graupel/rain collection efficiency; V G [= « l ( 7T ) 3 ] 1S mass ' 

weighted fall-speed for graupel. 


Pgacs = E as P'~\Vc ~ VsMsNog^HtA- + TfW + 

p f) A s Aq a S A G 


0.5 

A|Ag 


(A 27) 
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where Egs(~ 0.1) is the graupel/snow collection efficiency. 


„ naqc.EGcNoG ,Po^i B(b + 3) 

i GACW = : 1 — ; 2 T~Z 

4 p x b+3 


{A 28) 


where Egc(— 1) is the graupel/cloud water collection efficiency. 


Pw Acs = n>c.Esc a 


2 ^ nJi{ £»)i r(fc + (i) 


24 p p ' A^+ 6 


(A29) 


where n r (— is the number concentration of cloud water droplets; M c (— 4 x 10 12 ) is 

cv M c ' 

the average mass of cloud water droplet. 


Pgdep 


27rJVo G (Si-l) r 0.78 , ou a P\h/P Q \? r (^/ 2 + 5 / 2 ) i 
p{A" + B") [ X l V { (> A V2+5/2 J 

G 


(A30) 


, )T MT a -T)~ 1 

Pgfr = B^Nqr — ry i (A31) 

P X R 

where A^{= 0.66K -1 ) is the constant in Bigg freezing; B^(= 100m~ 3 s _1 ) is the constant 
in raindrop freezing equation. 
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Table Al List of microphysical processes and their parameterization schemes in Appendix 


Notation 

Description 

Scheme 

Pmltg 

Growth of vapor by evaporation of liquid water from the 
surface of graupel 

RH84 

Pmlts 

Growth of vapor by evaporation of melting snow 

RH83 

Prevp 

Growth of vapor by evaporation of raindrop 

RH83 

PlMLT 

Growth of cloud water by melting of cloud ice 

RH83 

PCND 

Growth of cloud water by the condensation of supersatu- 
rated vapor 

TSM 

Pgmlt 

Growth of raindrop by melting of graupel 

RH84 

PSMLT 

Growth of raindrop by melting of snow 

RH83 

Praci 

Growth of raindrop by the accretion of cloud ice 

RH84 

Pracw 

Growth of raindrop by the collection of cloud water 

R.H83 

Pracs 

Growth of raindrop by the accretion of snow 

RH84 

Praut 

Growth of raindrop by the autoconversion of cloud water 

LFO 

PlDW 

Growth of cloud ice by the deposition of cloud water 

KFLC 

PlACR 

Growth of cloud ice by the accretion of rain 

RH84 

Pi ho m 

Growth of cloud ice by the homogeneous freezing of cloud 
water 


Pdep 

Growth of cloud ice by the deposition of supersaturated 
vapor 

TSM 

PsAUT 

Growth of snow by the conversion of cloud ice 

RH83 

PsACI 

Growth of snow by the collection of cloud ice 

RH83 

PSACW 

Growth of snow by the accretion of cloud water 

RH83 

PSFW 

Growth of snow by the deposition of cloud water 

KFLC 

PsFI 

Growth of snow by the riming of cloud ice 

KFLC 

PSACR 

Growth of snow by the accretion of raindrop 

LFO 

PsDEP 

Growth of snow by the deposition of vapor 

RH83 

Pgaci 

Growth of graupel by the collection of cloud ice 

RH84 

Pgacr 

Growth of graupel by the accretion of raindrop 

RH84 

Pgacs 

Growth of graupel by the accretion of snow 

RH84 

Pgacw 

Growth of graupel by the accretion of cloud water 

RH84 

Pw ACS 

Growth of graupel by the riming of snow 

RH84 

Pgdep 

Growth of graupel by the vapor deposition 

RH84 

Pgfr 

Growth of graupel by the freezing of raindrop 

LFO 
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Figure Captions 

Fig. 1 Time evolution of (a) vertical velocity (mb hour -1 ), and (b) zonal wind (ms -1 ) taken 
from the TOGA COARE for a 20-day period. Downward motion in (a) and westerly wind 
in (b) are shaded. 

Fig. 2 (a) Linear correlation and (b) RMS of the temperature (solid) and specific humid- 
ity (dashed) between observation and simulation. The units of RMS in (b) are °C for 
temperature and gkg -1 for specific humidity respectively. 

Fig. 3 (a) Components of mass-weighted mean heat budget versus local change of mass- 
weighted mean temperature ( °Cday -1 ) and (b) Components of column-integrated mois- 
ture budget versus local change of column-integrated specific humidity (precipitable water) 
(mmday -1 ). Closed dots denote sum of mass- weighted mean latent heat of condensation 
and vertical potential temperature advection in (a) and sum of vertically-integrated con- 
densation (precipitation) and vertical moisture advection in (b). Open dots represent 
imposed large-scale horizontal temperature advection in (a) and moisture advection in 
(b). Symbols ”x” are surface sensible heat flux in (a) and surface evaporation flux in (b). 
Symbols delta denote radiative processes. 

Fig. 4 (a) Mass-weighted mean latent heat of condensation versus vertical potential tem- 
perature advection ( °Cday -1 ), and (b) vertically-integrated condensation (precipitation) 
versus vertical moisture advection (mmday -1 ). 

Fig. 5 (a) Sum of mass-weighted mean latent heat of condensation and vertical potential 
temperature advection (°Cday -1 ) versus surface rain rate (mmday -1 ), and (b) sum of 
vertically-integrated condensation and vertical moisture advection (mmday -1 ) versus rain 
rate (mmday -1 ). 

Fig. 6 As Fig. 3a except for those at (a) 194 mb, (b) 309 mb, (c) 487 mb, (d) 694 mb, and 
(e) 907 mb, and symbols ”x” denote vertical heat flux convergence. 

Fig. 7 As Fig. 3b except for those at (a) 194 mb, (b) 309 mb, (c) 487 mb, (d) 694 mb, and 
(e) 907 mb, and symbols ”x” denote vertical moisture flux convergence. Note different 
plotting scales in (a) and (b) . 

Fig. 8 (a) Column-integrated cloud content versus rain rate in (a)-(c), and (d) vertical profiles 
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of time- mean cloud contents. 


Fig. 9 (a) Column-integrated cloud microphysical processes related to conversions between 
vapor and cloud contents versus rain rate, and (b) vertical profiles of time-mean cloud 
microphysical processes. 

Fig. 10 Column-integrated cloud microphysical processes related to conversions between 
cloud water and other cloud contents versus rain rate. 

Fig. 11 Column- integrated cloud microphysical processes related to conversions between 
raindrop and other cloud contents versus rain rate. 

Fig. 12 Column-integrated cloud microphysical processes related to conversions between 
cloud ice and other cloud contents versus rain rate. 

Fig. 13 Column-integrated cloud microphysical processes related to conversions between 
snow and other cloud contents versus rain rate. 

Fig. 14 Column-integrated cloud microphysical processes related to conversions between 
graupel and other cloud contents versus rain rate. 

Fig. 15 Vertical profiles of time-mean cloud microphysical processes of (a) cloud water, (b) 
raindrop, (c) cloud ice, (d) snow, and (e) graupel. 

Fig. 16 Mean cloud microphysics budget. Units for cloud contents and conversions are mm 
and mm day -1 respectively. 
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